Draft version February 5, 2008 

Preprint typeset using 1^1^^ style cmulatcapj v. ^jllj^A, 



THREE-BODY KICK TO A BRIGHT QUASAR OUT OF ITS GALAXY DURING A MERGER 

LoREN Hoffman^ and Abraham Loeb^ 

Draft version February 5, 2008 

ABSTRACT 

The quasar HE0450-2958 was recently discovered (jMaeain at alJ2005|) to reside ^ 7 kpc away from a 
galaxy that was likely disturbed by a recent merger. The lack of a massive spheroid of stars around the 
quasar raised the unlikely suggestion that it may have formed in a dark galaxy. Here we explain this 
discovery as a natural consequence of a dynamical kick imparted to the quasar as it interacted with a 
binary black hole system during a galaxy merger event. The typical binary stalling radius provides a 
kick of order the escape velocity of the stellar spheroid, bringing the quasar out to around the observed 
radius before it turns around. This is consistent with the observed low relative velocity between the 
quasar and the merger-remnant galaxy. The gas carried with the black hole throughout the three-body 
interaction fuels the quasar for the duration of its journey, '--^ 2 x 10'' years. Gravitational radiation 
recoil could not have produced the required kick. 

Subject headings: black hole physics - galaxies: evolution - galaxies: interactions - galaxies: nuclei — 
(galaxies:) quasars: individual (HE0450-2958) - galaxies: starburst 



1. INTRODUCTION 

Observations indicate that every massive galaxy con- 
tains a supermassive black hole in its nucleus, which 
shines as a quasar during episodes of rapid mass ac- 
cretion. The vigorous gas inflows needed to fuel 
the most luminous quasars are thought to be trig- 
gered by galaxy mergers, in which tidal interac- 
tions channel large amounts of gas into the cen- 
tral ~ 100 pc fBvrd et al. 1987; Hernquist 1989). 
The maximum luminosities of quasars are strongly 
correlated with the properties of the ir host galaxies 
l|Shields et all l200,'^ IPTovd et alJ 120041) . in accordance 
with the tight correlations observed between black hole 
mass and stellar bulge mass (Magorrian, ct al. 199Sj 
iMarco ni fc Hunt 2003: Peng et al. 2005) and velocity dis- 
nersion (.Ferrare^^ fc Merritt 2000: Gebhardt et a,ljl200?i 
iTremaine et al.ll2002|) in quiescent galaxies. 

Given thes e correlations, the recent discovery 

l|Magain et al.ll2005D of a bright {My = -25.8) quasar 
at redshift z = 0.285 apparently not surrounded by a 
massive host galaxy is surprising indeed. For emission at 
the Eddington luminosity, the inferred black hole mass 
is ~ 4 X 10® Mq. There is a neighboring, luminous star- 
burst galaxy at a separation of ^ 1.7" corresponding 
to 7.2 kpc for a flat universe with a present-day Hub- 
ble parameter Hg = 71 km Mpc~^, and a matter 
density parameter = 0.27. This companion galaxy 
has a disturbed morphology and is heavily obscured 
by dust, typical of the Ultra-Luminous Infrared Galax- 
ies (ULIRGs) systematically associated with near equal- 
mass mergers (.Sanders et al. 1988 ) . Its visual magni- 
tude is My = —23. Close to the quasar lies a bright 
blob about 2 kpc across, consisting of gas ionized by 
the quasar radiation. The spectrum of the blob shows 
no stellar component and no dust extinction; hence the 
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host galaxy cannot be hidden by dust. The typical host 
for an My = —25.8 quasar is an elliptical galaxy with 
My « —23.2, about the magnitude of the neighboring 
ULIRG. However after carefully deconvolving a new HST 
image of the system, Magain et al. find no sign of any 
such galaxy surrounding this quasar. Fixing the half- 
light radius at — 10 kpc (a typical size for the host 
of an My ~ — 26 quasar), they put a conservative up- 
per limit of My = —21.2, 3.7cr fainter than expected, 
on the magnitude of any host. Alternatively, adopting 
Afy — —23 for the host magnitude, they set an upper 
limit of ~ 100 pc on R^, much too small for a normal 
galaxy of this magnitude. 

We suggest that the quasar observed in HE0450-2958 
may be a black hole that was ejected from the nucleus of 
the companion galaxy during the major merger causing 
its disturbed appearance. The gas now fueling the quasar 
was dragged with it from the galactic center when it was 
kicked out. If one of the merging galaxies harbored a 
binary black hole, then when the black hole from the 
other galaxy entered the system a complex 3-body inter- 
action would ensue, shortly ending with the ejection of 
the lightest bla ck hole a t a velocity of order the binary 
orbital velocity l)Heggiell 975^ . This scenario implies that 
one of the merging galaxies was relatively gas-poor (per- 
haps a giant elliptical), to allow survival of the binary 
against gasdynamical friction for a Hubble time before 
the merger. On the other hand at least one of the galax- 
ies must have been gas-rich, in order to fuel the quasar. 
The two galaxies must have been close in mass, to ex- 
plain the ULIRG appearance of the merger remnant and 
the high inferred black hole mass of the quasar, since it 
is the lightest body that gets ejected. 

The Eddington accretion rate for a black hole of mass 
M = 4 X 10** Mq is MEdd = (47rGMmp)/(ecrTc) = 
8.9eZjMQ yr~^, where rUp is the proton mass, ctt the 
Thomson cross-section, and e = 0.1e_i is the radiative 
efficiency. It takes ~ 2 x 10^ yr for the black hole to 
reach a turnaround radius of 7.2 kpc, so in order to re- 
main shining with e_i ~ 1 at the observed distance the 



2 



black hole must drag at least ~ 36% of its own mass 
with it in gas when it is ejected. The distances of clos- 
est approach in our three-body interactions are often 
larger than the radius conta ining this mass in an a-disk 
ijShapiro fc Tcukolskv 1983), and so the inner accretion 
disk can follow the black hole throughout the interaction 
and get ejected with it. 

2. BLACK HOLE BINARIES AND EJECTIONS 

When two galaxies of comparable mass merge, their 
black holes sink to the center by dynamical friction 
and form a gravitationally bound binary. Once the bi- 
nary separation reaches the "hardening radius", rhard ~ 
G[min{Mi,M2)]/4:a^ ~(l-5) pc for massive galaxies 
(where a is the ID stellar velocity dispersion and Afi.2 
are the black hole masses), its further reduction is dom- 
inated by close encounters with stars that pass within 
the binary orbit. These stars undergo strong 3-body in- 
teractions with the binary and escape its vicinity with 
speeds comparable to the binary orbital speed. Only 
stars on nearly radial orbits can contribute to the hard- 
ening of the binary in this stage. In the low-density nuclei 
of massive elliptical galaxies, the total mass in stars on 
these so-called "loss cone" orbits is small compared to 
the mass of the binary. Since the two-body relaxation 
time is long compared to the Hubble time, stars cannot 
diffuse back into the loss cone as fast as they are kicked 
out. If all stars on loss cone orbits run out, the decay 
of the binary ceases. Since the binary ejects of order its 
own mass in stars per e- folding in radius, one naively ex- 
pects binaries in the cores of massive elliptical galaxies 
to stall at separations of order 1 pc. If the separation 
reaches ~0.01-0.1pc, then the binary would quickly co- 
alesce through the emission of gravitational radiation. 
The question of whether and how it crosses the gap from 
'--^ 1 pc to this smaller scale is known as t he "final parsec 
problem" ( M erritt fc Milosavlievi3l2004j) . 

Black hole binaries may coalesce in some cases but 
not others. For example, the mass of any close binary 
companion to the nuclear black hole of the Milky Way 
is limited to M2 S 10'* by the Keplerian orbits of 
stars near the galactic center (jMerritt fc Milosavlievid 
l2004jl . The so-called X-shaped radio galaxies have been 
explained as jets bent by sudden changes in the spin of 
a black hole vipon coalescence with a binary companion 
(^crritt & Milosavlicvic 200^. Gas drag may lead to 
coale scence in gas-rich mergers (jMerritt fc Milosavlievid 
l2004j) : however the gas-poor cores of massive ellipticals 
could allow binaries to stall. There are some tentative 
observational hints that sub-parsec scale binaries exist, 
including periodic variations and double-peaked broad 
line profiles in active galaxies, and central minima in the 
surface brightness profiles of some early type galaxies 
which are possibly associated wi th the clearing of the 
loss cone around stalled binaries ('Komossa"2003^ . 

When an unequal mass black hole binary coalesces by 
gravitational radiation, the waves carry off a net linear 
momentum and t he merged bl ack hole recoils in the op- 
posite direction l)Favata et al J 12004'). However 300 km 
s""'^ is a reasonable upper li mit to t he velocity achievable 
by th is mechanism ( Merrit t" et al.l 1^004: Blanchet ct a3 
mass ejected radially at 300 km s ^ from the 
center of the average stellar spheroid hosting a 4 x IO^Mq 
black hole (M^,„ = 1.7 x lO^M©, Re = 3.6 kpc) would 
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Fig. 1. — Black hole trajectories following a 3-body ejection. The 
blue (thick) line shows the path of the intruder and the red and 
yellow (thin) lines that of the binary members; each is plotted only 
until its first return to the galactic center. Time labels in units 
of 10^ yrs are marked off at various points along the trajectory, 
going out and returning. In the upper left inset we zoom in on the 
innermost ±7 pc. Here the trajectories are plotted only up to the 
point of ejection (with return trajectories left out). One can clearly 
see the binary and intruder flying off in opposite directions at the 
end of the interaction. The lower panel shows the distribution of 
ejection velocities over all ~10^ runs. 



only reach a distance of ~ 150 pc before turning around, 
far smaller than the separation of the quasar in HE0450- 
2958 from its companion galaxy. 

If a binary stalls and its host galaxy merges with an- 
other galaxy, a third black hole may spiral in and undergo 
strong three-body interactions with the pair. The merger 
event will drive mass through the binary at a significant 
rate, partially refilling the loss cone, but the final stages 
of the merger transpire so rapidly that the intruder is 
likely to arrive in the vicinity of the bi nary before the 
binar y has a chance to coalesce (Milosavli evic fc Merritt! 
1200 1"). The three-body interaction will end in the harden- 
ing of the binary by an amount comparable to its initial 
binding energy, and the ejection of the lightest bo dy at a 
speed comparable to the binary's orbital speed l)Heefgid 
119751) . The binary's orbital speed at the stalling radius is 
a few times a, i.e. of order the escape velocity Vesc from 
the stellar bulge (e.g. for a Hernquist profile (Hernauis^ 
Il990|) Vesc = 2^/2veirc{r = fl) with a ^ Varc)- Thus, 
qualitatively we expect three-body kicks to yield ejection 
velocities of order the bulge escape velocity, resulting in 
many near-ejections. This is consistent with the large 
distance and low velocity of the HE0450-2958 quasar rel- 
ative to the ULIRG. Figure 1 shows an example of such 
a near-escape trajectory. 

3. NUMERICAL CALCULATIONS 

To make our dynamical arguments more quanti- 
tative, we ran a series of 83,250 3-body scatter- 
ing experiments using Svcrre Aarseth's KS-regul arized 
^Kustaanheimo & Stiefel 1965) Chain code tAar sethl 
120031(1 to compute the black hole trajectories. In each 
case the intruder's mass was set to M3 = 4 x 10® Mq and 
the masses of the initial binary members were chosen 
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uniformly in their logarithm between 1.0 and 3.5 times 
the mass of the intruder, a generic range for major merg- 
ers. From here on we will refer to the initial black hole 
binary as the inner binary and the "binary" formed by 
the intruder and the inner binary's center of mass as 
the outer binary. The semi-major axis abin of the in- 
ner binary was always initialized to r hard/ 4:, and that 
of the outer binary was set using th e stability criterio n 
of Mardhng and Aarseth (Mardh ng fc AarsethI IT999|) . 

(i?°"7a,„) « 2.8[{l + qout){l + eout)/VT^^h which 
reliably estimates the most distant intruder orbit for 
which strong 3-body interactions can begin. Here 
is the pericentre distance of the outer binary, a™ is the 
semimajor axis of the inner binary, qout = M3/{Mi-\-M2) 
is the mass ratio of the outer binary, and Cout is the outer 
binary's eccentricity. The initial eccentricity of the inner 
(outer) binary was chosen uniformly between 0.0 and and 
0.2 (0.3), where this range of low eccentricities naturally 
results from dynamical friction. Finally, the three Euler 
angles of the intruder's orbital plane were chosen ran- 
domly relative to the reference plane of the binary orbit, 
as was the phase of the initial periastron of the binary. 
Both orbits were always started at pericentre; since many 
orbital periods elapse before unstable interactions begin, 
the relative phase is effectively randomized in any case. 

The Chain code allows the user to add conserva- 
tive or velocity-dependent external perturbations, and 
so we added a dynamical friction force given by Chan- 
drasekhar's formula fCh andrasekhar 1943() to the mo- 
tion of each body. The density and Coulomb logarithm 
In A were chosen to give a dynamical friction timescale 
10^ yr. Each run was terminated when (a) any one 
of the bodies was expelled to a distance of 1 kpc from 
the center of mass; (b) any two bodies came within 30rs 
(where = 2GM/(? is the Schwarzschild radius) of each 
other or the gravitational radiation timescale dropped 
below 10^ yr (roughly the dynamical time) for any pair, 
or (c) the total time elapsed reached 6 x 10^ yr. We 
selected out just those runs where the intruder never 
came within 5000 Schwarzschild radii (~0.2 pc) of ei- 
ther bi nary member , reflectin g the scale below which an 
a-disk ijShapiro fc Tcukolskvl ll98^ with sufficient mass 
to fuel the quasar might be disrupted. This is also 
about the size inferre d for the broa d-linc region of the 
HE0450-2958 quasar ijCreene fc Holl2005t iMerritT'et aO 
linnS). Of the runs not terminated by criterion (b) or 
(c), the closest approach of the intruder was > 5 x lO'^r^ 
in 36% and > lO^'r^ in 22% of the runs. To each 
successful run we then assigned a galaxy based on the 
empirical black hole/bulge correlatio ns. Each galaxy 
was modeled as a Her nquist profile iHernauisti Il990() 
with a bulge virial mass (jM arconi fc HunHl2nn3D « 
3Re(7^/G chosen according to a log-normal distribution 
with a mean given by the relation (Marconi fc Hunt 
l2t)03l) ~ (1.71 X 10"Mr;,)fM/( 4 X 10^ M^)]^-"-^ 

and standard deviation of 0.25 dex ijMarconi fc H^mtl 
llOOS). The v elocity dispersion was calculated from the 
M-cr relation ijTremaine et al.ll2002|) . a « (262 km s'^) 
[M/(4 X 10^Mo)]°-249, and the bulge scale radius set by 
Re = GMmr/3a^. We then added to the stellar bulge a 
dark matter halo with a mass ratio that is an order of 
magnitude larger than the cosmic matter-to-baryon ra- 
tio (corresponding to a star formation efficiency of 10%) 
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Fig. 2. — Probability P of observing the black hole at various 
projected separations D and line-of-sight velocities v. The colors 
are scaled in log P(D,?;), where P(D,v) is the probability density 
for obtaining the values {D,v). The probability incorporates the 
length of time spent by the ejected black hole at these values, av- 
eraged over ~ 10'' successful runs (which were neither terminated 
before the ejection nor eliminated by our minimum-distance cri- 
terion). The white star shows the estimated coordinates of the 
HE0450-2958 quasar based on the radial velocity range reported 
in Magain et al. , with errors bracketed by the white lines. The 
yellow (thicker) lines indicate the magnitude of the quasar's red- 
shift relative to the galaxy (30±20 km/s) measured by Merritt et 
al. (2005). The histograms show the corresponding projections of 
P onto the D and v axes. 



and a concentration parameter of 9.Q( Mhnio I 2.10 x 
10^Mq)-°13 with a scatter of 0.18 dex l|Bullock et alJ 

l2ml . 

The black hole was ejected radially from the galac- 
tic center with the final velocity from the 3-body run 
and its trajectory integrated for 6 x 10^ yr. We im- 
posed dynamical friction with In A = 3, which had al- 
most no effect on the trajectories in high velocity ejec- 
tions, but efficiently damped out those with low ve- 
locities. Finally, we chose a random line of sight and 
projected all distances and velocities accordingly. We 
added up the total time spent in each logarithmic dis- 
tance and velocity bin to assess the relative probability 
of seeing the quasar at different projected separations 
and line-of-sight velocities. Figure 2 shows the result- 
ing probability distributions. The observed parameters 
of the quasar under consideration, D ~ 7.2 kpc and 
V < 200 kms~^, appear to be fairly typical in the statis- 
tics of our runs. Note that the results scale simply with 
the black hole mass in our model. If the inner binary 
always starts at a fixed fraction of the hardening radius, 

then Mbh cx a* gives abin oc Mbh/o''^ « . We then 

have Vej oc ^ Mbh/abin oc M^l^ cx a. The self-similarity 
of the Hernquist model implies that the turnaround dis- 
tance rmax in units of the galactic scale radius is fixed 
by Vej in units of a. Since Vej/u is constant, we get 

1/2 

fmax (x Re CK. M^^ . Hcncc neglecting minor deviations 
from self-similarity due to dynamical friction, we expect 

the velocities in our plots to scale roughly as M^j^'^ and 

distances as M^^^ . 

4. DISCUSSION AND CONCLUSIONS 

We have proposed that the HE0450-2958 quasar, 
around which no galactic host has so far been detected. 
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was ejected from the nearby galaxy in a three-body in- 
teraction following a merger of a massive gas-rich galaxy 
with a giant elliptical. This scenario explains (a) the 
fueling of the bright quasar displaced by kpc from 
the galactic center; (b) the survival of the binary in the 
gas-poor galaxy before the merger; and (c) the disturbed 
appearance of the companion galaxy. One may also ask 
how the presence of the "blob" and narrow line emis- 
sion in the quasar spectrum, both indicative of ionized 
gas in the quasar's vicinity far outside the black hole 
radius of influence, fit into the ejection picture. Mer- 
ritt et al. (2005) infer a size of roughly 1.5 kpc for the 
narrow- line region (NLR). Since the gas at this distance 
cannot possibly be bound to the black hole, they argue 
that the very small observed redshifts between quasar, 
galaxy, and narrow emission line gas rule out the ejec- 
tion scenario. However since the ejected quasar would 
spend most of its time near turnaround, we find a quite 
small velocity difference between the quasar and ULIRG 
to be consistent with the high-velocity ejection model. 
For instance, the black hole modelled in Figure 1, with 
Vej ^ 1200 km s~^, spends 17% of its time with v < 100 
km/s. After projecting onto the line of sight and aver- 
aging over many runs as in Figure 2, one can see that 
very small velocities are still prevalent in our statistics 
at distances around 7 kpc. 

In galaxy merger simulati ons including feedba ck from 
star formation and AGN l|Springel et alJ l2005|) . large 
amounts of gas are blown from the galactic center out 
to distances of ^5-10 kpc during the coalescence of the 
two nuclei. The disturbed system also shines as a ULIRG 
during most of this phase. Once the feedback shuts off, 
the gas that does not escape settles back into the in- 
ner galaxy on the dynamical timescale of ~10® yrs. The 
mass in narrow line emitting gas inferred from absolute 
line fl uxes in typical NLR is around 1O^~^M0 (Peterson 
^^3)- Ten million years after the core merger, more 
than this mass in gas below 10^ K would likely be avail- 
able within the few kpc sphere surrounding a black hole 
ejected to a distance of ^7 kpc. It is plausible that 
clouds formed from this cool gas could produce the ob- 
served narrow line emission in the quasar spectrum while 
the black hole (and gas) are near turnaround. The ob- 
served blob would be a portion of the gas photoionized 
by anisotropic radiation from the quasar. It is not feasi- 
ble for Bondi accretion from this external gas to fuel the 
quasar, as the densities required are unrealistically high. 
Rather the quasar must be fueled by the remnant inner 
accretion disk dragged with it from the galactic center, 
which can remain bound to it through the ejection as 
demonstrated above. 

Our proposed scenario is viable as long as the merger 
process does not induce coalescence of the binary through 
dynamical friction on a rearranged distribution of stars 
and gas, before the three black holes interact. For typical 



parameters, we find that the feeding of the loss cone of 
stars by tidal interactions during the merger is not suffi- 
cient to cause coalescence of the binary before the third 
black hole sinks in (Roos 1981). Coalescence through 
friction on gas would be possible by momentum conser- 
vation only if a large quantity of gas, comparable to or 
larger than the binary mass, is channeled through the 
interior of the binary before the three black holes inter- 
act. However, since the gas initially resides within the 
quasar's original host galaxy, it is more likely to fuel the 
quasar before reaching the nucleus of the dry elliptical 
that hosts the binary. Feedback from the quasar could 
then easily supply sufficient energy to prevent the gas 
which gets t idally-stripped from th e outer envelope of its 
host galaxy l|Wvithe fc Loell 20031 from settling into the 
few parsec scale region surrounding the binary before the 
three black holes interact. 

Both our calculations and the identification of H E0450- 
2958 as a quasar without a typical host iMagai net alJ 
120051) are based on the correlations between black hole 
masses and their host bulges established for isolated, 
quiescent galaxies. The extension of these correlations 
to active and me rging system s is somewhat less certain. 
Severa l authors ijGreene fc H o 2005; Onken et al. 200| 
[Wandcl 2002) have studied the Mbh-o- and MBH-Mbuige 
relations in active galaxies and reported statistically sig- 
nificant deviations from the slope and zero-point mea- 
sured for inactive systems, as well as slightly larger scat- 
ter. However the reported deviations are not typically 
greater than a factor of 2-3 and hence qualitatively sup- 
port our application of the correlations to an interacting 
system. 

Since the mass of the HE0450-2958 black hole and cor- 
responding expected host properties are uncertain, one 
cannot prematurely conclude that it is truly a "naked" 
quasar. Further observations, particularly in bands ad- 
mitting a higher ratio of galaxy to quasar luminosity, are 
needed to place tighter constraints on the properties of 
any host galaxy. This paper demonstrates the plausi- 
bility of a three-body ejection producing a naked quasar 
near a recently-merged galaxy. If such a system is unam- 
biguously observed, this would provide a key link in our 
understanding of galaxy formation and the hierarchical 
buildup of supermassive black holes. 
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